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Abstract 

This paper considers the Cramer-Rao lower Bound (CRB) for the source locaHzation problem in the near field. 
More specifically, we use the exact expression of the delay parameter for the CRB derivation and show how this 
'exact CRB' can be significantly different from the one given in the literature and based on an approximate time 
delay expression (usually considered in the Fresnel region). This CRB derivation is then generalized by considering 
the exact expression of the received power profile (i.e., variable gain case) which, to our best knowledge, has 
been ignored in the literature. Finally, we exploit the CRB expression to introduce the new concept of Near Field 
Localization (NFL) region for a target localization performance associated to the application at hand. We illustrate 
the usefulness of the proposed CRB derivation and its developments as well as the NFL region concept through 
numerical simulations in different scenarios. 
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I. Introduction 

Sources localization problem has been extensively studied in the literature but most of the research works are 
dedicated to the far field case, e.g., ifTl. |[2l. 

In this paper, we focus on the situation where the sources are in a near field region which occurs when the source 
ranges to the array are not 'sufficiently large' compared with the aperture of the array system Q, ||4l. Indeed, 
this particular case has several practical applications including speaker localization and robot navigation [51, |[6l . 
underwater source localization Q, near field antenna measurements lU, 121 and certain biomedical applications, 
e.g., ifTOl . Recently, some works considered both far filed and near filed localization ifTTI . |[T2l where the authors 
propose different methods to achieve a better localization performance when the source moves from far field to 
near field and vice versa. 

More specifically, this paper is dedicated to the derivation of the Cramer Rao Bound expressions for different 
signal models and their use for the better understanding of this particular localization problem. CRB derivation 
for the near field case has already been considered in the literature [Si, |[T3l . |[T4l . |[T5l . In Q, |[T3l . the exact 
expression of the time delay has been used to derive the unconditional CRB in matrix form, i.e., expressed and 
computed numerically as the inverse of the Fisher Information Matrix (FIM). In |[T4l . the conditional CRB based 
on an approximate model (i.e., approximate time delay expression as shown in Section ITl-AI ) is provided. 

Recently, EL Korso et al. derived analytical expressions of the conditional and unconditional CRB for near field 
localization based on an approximate model |[T5l . In the latter work both conditional and unconditional CRB of 
the angle parameter are found independent from the range value and are equal to those of far field region. 

In our work, we propose first, to use the exact time delay expression for the derivation of the conditional and 

unconditional CRB and provide closed form formula that are compared to those given in [151. The development 

(i.e., Taylor expansion) of the exact CRB allows us to highlight many interesting features including: (i) a more 

accurate approximate CRB for the exact model as compared to the CRB based on an approximate model, (ii) and 

a detailed analysis of the source range parameter effect on its angle estimation performance , 

Secondly, we take into consideration the spherical form of the wavefront into the power profile. Indeed, in the 
'Part of this work has been published in 1161 



near field case, the received power is variable from sensor to sensor which should be taken into account in the 
data model. By considering such variable gain model, we generalize the previous CRB analysis and investigate the 
impact of the gain variation onto the localization performance limit. 

Finally, we propose to exploit the exact CRB expression to specify the 'near field localization region' based on 
a desired localization performance. In that case, the 'near field localization region' is shown to depend not only on 
the source range parameter and array aperture but also on the sources SNR and observation sample size. 

The paper is organized as follows: Section |ll] introduces the data model and formulates the main paper objectives. 
Section [nil the exact conditional and unconditional CRB derivations and their Taylor expansions are provided for 
the equal gain case. Section |IV] generalizes the previous analysis to the variable gain model and highlights the 
differences between equal and variable gain cases. In Section |Vl we introduce the concept of near field localization 
region and illustrate its usefulness through specific examples. Section |Vl] is dedicated to simulation experiments 
while Section IVIII is for the concluding remarks. 

II. Problem formulation 

A. Data model 

In this paper, we consider a uniform linear array with N sensors receiving a signal, emitted from one source 
located in the near field region, and corrupted by circular white Gaussian noise w„ of covariance matrix cj^Iat. The 
^th ^]-j-^y output, n = 0, ■ ■ ■ , — 1, is expressed as 

Xn{t) = s{t)e^^- +Vn{t) t = l,---,T (1) 

where s{t) is the emitted signal. The exact expression of the time delajo r„ is given by 



27rr / / n'^cP 2nds'm6 



r. = — h/l + -3 ; 1 (2) 



where d is the inter-element spacing, A is the propagation wavelength and (r, 9) are the polar coordinates of the 
source as shown in Fig. [1] In the literature, the near field region (also called Fresnel region) is given by lITTl 

0.62 ( ^-j^ j < r < 2^^^^ (3) 



The first sensor, n = 0, is considered for the time reference. 
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Fig. 1. Near field source model 

Most existing works on near field source localization consider the following approximation of the time delay 
expression to derive simple localization algorithms as well as CRB expressions, e.g., ifTSl |[T5l 



In the sequel, the source signal will be treated either as deterministic (conditional model) or stochastic (uncon- 
ditional model). Indeed, in the array processing, both models can be found: 

1) Conditional model in which the source signal is assumed deterministic but its parameters are unknown. 

2) Unconditional model in which we assume that the source signal is random. In our case, we will assume s(t) 
to be a complex circular Gaussian process with zero mean and unknown variance a^. 



1) Based on the exact expression of the time delay Tn in equation Q, we aim to derive the exact conditional and 
unconditional CRB and compare them with the existing CRB given in ifTSl and derived from the approximate 
model in (4). 

2) The distance from the source to the n*'* sensor dn is a function of the sensor position (i.e., sensor index) 
according to 




(4) 



B. Objectives 



d. 



'n 



\J r"^ — 2ndr sm{9) + n^d^ 




(5) 
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where r is the distance from the source to the reference sensor (i.e., for n = 0, do = r). Hence, the 
received power profile is variable from sensor to sensor. In the far field case, this variation is negligible but 
not necessarily in the near field context. To the best of our knowledge, this point has not been taken into 
consideration in the existing literature. We would like to investigate the impact of such gain variation into 
the localization performance limit. 
3) The near field region has been so far assimilated to the Fresnel region which depends on the antenna size 
and the signal wavelength only. However, the localization performance depends on other system parameters 
(SNR, sample size, angle position, • • • ) for which reason we introduce the concept of near field localization 
region where the localization error is upper bounded by a desired threshold value depending on the considered 
application. 

HI. Conditional and unconditional CRB derivation with Equal Gain (EG) 

Here, we do not consider the variation of the power profile which is equivalent to assuming that all sensors 
have equal gain as shown in equation (1). The goal of this Section is to derive the conditional and unconditional 
CRB for the angle and range parameter estimation of near field source and to make comparison with the results 
presented in |[T5l . 

A. Exact conditional CRB with EG 

In this Section, we consider the source signal as deterministic according to the model 

sit) = a(t)e^'(2-/o*+'^W) 

where /o is the known carrier frequency while a{t) and tlj{t) are the unknown amphtude and phase parameters of 
the source signal. Under the data model assumption of Section III-AI we derive next the exact deterministic (i.e., 
conditional) CRB for the location source parameter estimation. 

1 ) Exact conditional CRB derivation: In this deterministic case, the log-likelihood function of the observations 
is given by 

= -NTlrnr - NTlna^ - \\x - /xlP 



where 



X = 


[x^(i)... ,x^(r)f 


x^(i) = 


[xo{t), ■ ■ ■ ,XN-l{t)f 


n = 


[s(l)a^(0,r),... ,s{T)a^{e,r)f 


a{e,r) = 




i = 




* = 




Q = 


[«(!),••• ,aiT)f 



where ||.{| refers to the Frobenius norm and [.]-^ is the transpose operator. 

The CRB is equal to the inverse of the Fisher Information Matrix (FIM) defined by 



the latter is given in our particular case by 



[FIM(^)],,,- = — — — + —TZe{-^^} 



(6) 



(7) 



where TZe{.} refers to the real part of a complex valued entity. The FIM is given by the block diagonal matrix as 

Q 



FIM 



27V T 



NT 



Indeed, it is shown in |[T5l . that the FIM of the desired location parameters is decoupled from the noise variance 
cr^ and the sources magnitude parameters ex. Consequently, the CRB of the range and angle parameters is equal to 
the 2x2 top left sub-matrix of the inverse matrix 
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frd — fer — 2TD3jyf]^(T g fr) 
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89''"' 86 



Moreover, the vectors of size T x 1, f^g, fL,, f^^ and fL, are given by 



where represents the Hadamard product and Iat is the all-one vector of size N x 1. Finally, matrix F^^ of size 
T X T, is given by 

2iV , , 
= — Trdiagia. ot) 

which is the diagonal matrix formed by vector a Q a. 

By using the block structure of the FEVI and the Schur's matrix inversion lemma |[T9l . we obtain the closed form 
expression of the CRB given by the following lemma: 

Lemma 1: The non-matrix expressions of the exact conditional CRB in the equal gain case for a source in the 
near field, for > 3 and 9 7 ^ ifc^, are given by 



CRB^,(r) = 
CRB^,(0) = 







2™snr) 


E,,{e)E,,{r)~E,,ir,er 






2™snr) 


E^,{e)E,,{r)-E^,{rfiY 


(2^'^SNR 






) E,,{e)E,,{r)-E^,(r,er 



(8) 
(9) 
(10) 



where CRBg (r, 6) is the non diagonal entry of the considered 2x2 CRB matrix (it represents the coupling between 
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the 2 parameters) and 

Eeg{e) = Wref - ^{llrgf 



Proof: The proof of this lemma can be deduced directly from lemma |5] proof given in appendix lAl 

2) Taylor Expansion (TE) of the CRB: This section aims to highlight the effects of some system parameters on 
the localization performance and to better compare our exact CRB with the CRB given in ifTSl . For that, we propose 
to use a Taylor expansion of the expressions of lemma [T] For simplicity, we omit the details of the cumbersome 
(but straightforward) derivations and present only the final results in the following lemma: 
Lemma 2: The Taylor expansions of the exact CRB expressions of lemma [T] lead to 

3A2 



V2{N)-QiN -iMioN"^ -lf>N +ll)sm{0)- 

r 

1 (P 
+—(2N - l)(384iV3 - 1353iV2 + 1379iV - 368)^ 
70 

1 

+— (186A^^ - 1590iV3 + 53517V2 - 6795iV + 2890) W{e)-;^ 



(11) 



CRB^,(r) 



6r^A 



2 \2 



TD 



SNR 



(i47r2 COs4(e)p3(iV) 



[I5r2 - m{N - l)s\n{e)dr 



+— {sin2(6l)(1061iV2 - 2625iV + 2911) 



+225iV2 - 315iV - 135} (f 



(12) 



where p2{N) = {8N - ll)(2iV - 1) and p^iN) = N{N^ - 1){N'^ - 4). 

These expressions are useful to compare our CRB to the one in |[T5l but also, they allow us to reveal how the 
different system parameters affect the localization performance. For example, one can see in particular that: 

• The CRB decreases linearly w.r.t. the observation time as well as the deterministic SNR. 

• The CRB goes to infinity when 9 ^ ±^ and the best localization results are obtained in the central direction, 
i.e., for (9 = 0. 



• Asymptotically (for large antenna sizes), the first term of the TE decreases in N'^, i.e., if we double the number 
of sensors, the CRB will be decreases approximately by a factor of 8. 

3) Comparison with the 'CRB of the approximate model': Let us first recall the CRB expressions given in |[T5l 
and based on the approximate model in (01): 

Lemma 3: The non-matrix expressions of the approximate conditional CRB |[T5l for a source in the near field, 
for > 3 and 9 ^ ib|, are given by 

= ^9f^ -ir-2 27^wM^2(iV) (13) 

PRO ^ 6^^^^ 

(I5r2 + mdr(N - 1) sin(0) + dV(A^) sin2(0)) (14) 



Comparing lemma |2] to lemma |3l one can make the following observations: 

• First, we note that the first and main term of equation (fTTI ) (resp. of equation (fT2l )) is equal to the first term 
of equation ([T3] ) (resp. of equation (O). 

• We note also, that the Taylor expansion (of the time delay) followed by CRB derivation leads to different 
results as compared to CRB derivation followed by Taylor expansion (of the CRB). The latter expansion being 
more accurate than the former as illustrated by our simulation results. 

• The approxim.ate CRB of the angle estimate shown in lemma |3] is independent of the range parameter (it is 
the same as the far field CRB) while the expression of CRB(^) in lemma |2] reveals how it is affected by the 
range parameter. In particular, at the first order, one can see that the CRB decreases (resp. increases) as a 
function of ^ for 6* € [0, f [ (resp. for 6^ g] - f , 0]). 

B. Exact unconditional CRB with EG 

For unconditional CRB, the signal is assumed Gaussian complex circular with zero mean and variance a^. 
Assuming that ^ = [0, r, cjg , cr^]^ is the vector of the unknown parameters, the log-likelihood function of the 
observed data is given by 

ln(|) = -A^Tlnvr - In detCE) - Ttr(5]-iR) 
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where "E = a1a{6,r)a{9,r)^ + cj^Iat (theoretical covariance matrix), R = ^ Y^=i x{t)x^ (t) (sample estimate 
covariance matrix), tr(.) refers to the matrix trace and (.)^ is the transpose conjugate operator. 
The Fisher Information matrix in this Gaussian case is given by 



[p,m«„.„^t..(e-.||e-.|) 



(15) 



Since the CRB is equal to the inverse of the Fisher information matrix and using the results of 11201 for the two 
desired localization parameters, we have the following lemma: 

Lemma 4: The non-matrix expressions of the unconditional CRB in the equal gain case for a source in the near 
field, for > 3 and 6* / ±^ are given by 

:i + (SNR) iV)I)sNR" 



CRB^f^C^) 
CRB:f,(r) 



(SNR)2 N 
:i + (SNR) iV)I)sNR' 



(SNR)2 N 
(1 + (SNR) iV)i?SNR' 



CRB^3(r) 
CRB^,(r,^ 



(16) 
(17) 
(18) 



(SNR)2 N 

where SNR = ^ and £'snr x CRBg^ represents the normalized CRB^^ (see lemma[T]) depending on the localization 
parameters, the array geometry, and the sample size only. 

Proof: The proof of this lemma can be deduced directly from lemma |6] proof given in appendix (iBl ). 

This result translates the fact that the unconditional CRB varies in a similar way as the conditional CRB w.r.t. 
variables 9, r and T. 

However, concerning the SNR parameter and the number of sensors, it is interesting to observe that: 

• At low SNRs (i.e., if SNR x N << 1), the CRB decreases quadratically (instead of linearly in the conditional 
case) w.r.t. the SNR and in (instead of N^) in terms of the number of sensors. 

• However, for large SNRs (i.e., if SNR x N » 1), the unconditional CRB behaves similarly as the conditional 
CRB w.r.t. parameters SNR and A'^ (i.e., it decreases linearly with the SNR and in w.r.t. the number of 
sensors). 
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IV. Conditional and unconditional CRB derivation with Variable Gain (VG) 

In this case, we consider the same model as that described in Section III-AI except that the received power is 
variable from sensor to sensor. The n^^ output array is expressed as: 



xn{t) = 7n(r, e)s{t)e^^-^^''^ + vn{t) = s(t)(7n(r, 0)e^'^"(^'^)) + Vn{t) t = l,. 

where 7n(?', 0) represents the power profile of the n^^ sensor and is given by ifTTl 

1 1 



(19) 



dn 



(20) 



dn being the distance between the source signal and the n*^ sensor. 

As we can see, in this model both the time delay profile and the power profile carry information on the desired 
source location (r, 6). Our objective here, is to investigate the roles of both profiles in the performance limit given 
by the CRB. 



A. Exact conditional CRB with VG 

Similarly to Section IIII-AI the conditional CRB is given by equation (7) which reduces in the variable gain case 



to 



FIM 



Q 







(21) 



which translates the fact that the FIM of the desired localization parameters is decoupled from the noise variance 
(T^ but not from the source magnitude parameter as in the equal gain case. Hence the CRB matrix of the range and 

angle parameters is equal to the 2x2 top left sub-matrix of the following inverse matrix 

-1 
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where fre = for and 



2™SNR(ll7ef + (7 if {re r,)) 



(22) 
(23) 
(24) 



frr = 2ri5sNR(||7,f + (7 7f (^rOr,)) 
/,e = 2r5sNR((7f7r) + (7 7f (t^Ot.)) 

Note that £'snr is defined here similarly to -D^nR ^^SNR ~ ^'T^^'^ implicitly scaled by the square of 

the range value which has been extracted from the received signal amplitude in ([T9l ) in order to explicit the role of 
the power profile. 

Consequently, ^^snR ^^^^ square meter (while i^SNR ^ constant without unit). Also, this way of normalizing 
the signals leads to CRB expressions in the variable gain case that are scaled (multiplied) by a factor of as 
compared to the equal gain caseQ For this reason, when comparing the equal gain scenario to the variable gain 
one in Section |Vll we normalize the received signal power by dividing ||q:||^ by in the former case. 

Moreover, the vectors of size T x 1, f^g, fj^, f^j. and f^^ are given by 





— ^9jp — 


2 

— {■yQjfre 




- - 

rip 


^(7 7)^Tr 


foe 


— *^9a — 


^(7^70)« 
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= F = 

ra 


^(7^7r)« 



where 



7 

70 

7r 



ho{r,0), ...,77v_i(r,6')]' 

89''"' 89 

8^0 8^N-i 



dr 



dr 



Finally, matrices F^^, Faa, ^ipa of size T x T, are given by 



a 



2 2 

"7II diag(Q;0Q;), F^^ 



2 „ „2 

II7II It and F^q, = 0. 



a 



This can be seen as if in the equal gain case, the received power is set equal to a constant value while in the variable gain case, it is the 
transmit power that is set equal to a constant. 
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Again by using the Schur's matrix inversion lemma, one can obtain: 

Lemma 5: The non-matrix expressions of the conditional CRB in the variable gain case for a source in the near 
field, for > 3 and 9 ^ are given by 

1 \ E,g{e) 



where 



CRBS,(r) 



CRBS,(r,( 



1 \ E^g{r) 



2TD 



SNR. 



2TD 



SNR/ ^^'9^ 



E,g{9)E,g{r) - E,g{r,ey 

Evg{r,d) 
E,g{e)E„g{r) - E,g{r,eY 



llTef + (7 7)^('re0re) 
llTrf + (7 7f (-TrOr,) 
7e7r + (7 7)^(Te0Tr) - 



[((7 7fr,)2 + (7^7.r 



1711 
1 



[((7 7fr.)2 + (7%.)^^ 



I7l 



7^ [(7 7)^-^0(7 ifrr + (7^7.) (7^70)] 
I7ll 



(25) 
(26) 
(27) 



(28) 
(29) 
(30) 



Proof: See appendix lAl 

One can observe, that the equal gain case (i.e., lemma [U is a particular situation of lemma |5] where 7 = l^v 
and ^SNR is divided by r^. 



B. Exact unconditional CRB with VG 

Under the stochastic Gaussian model of Section lTlI-Bi the unconditional CRB with variable gain can be expressed 
as stated by the following lemma: 

Lemma 6: The non-matrix expressions of the exact unconditional CRB in the variable case for a source in the 
near field, for > 3 and 6* 7^ are given by 



Proof: See appendix |B] 



CRBig[e) = 


( (i+SNR||7in^sNR^ 
I (SNR)^||7||^ - 


1 CRBigiO) 


(31) 


CRB^gir) = 


( (i+SNR||7in^sNR^ 
I (SNR)2||7||^ , 


1 CRBS,(r) 


(32) 






|CRBS,(r,0) 


(33) 


[ (SNR)^||7ir J 
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Similarly to the equal gain case, D^-^^ x CRBg^ represents a normalized CRB independent from the SNR value. 
Also, SNR is defined in the same way as the SNR in lemma |4] except that it has been implicitly scaled by the 
square of the range parameter (SNR unit is now square meter). 

C. Comparison between VG and EG cases 

By considering the extra-information given by the received power profile, one is able to achieve better localization 
performance as will be shown later in Section |Vll To better compare the CRB expressions in the constant and variable 
gain cases, we provide here the first term of their Taylor expansion 

3A2p2(iV) 



3AVp2(iV)(l + AViW) 



2^'™SNR^' cos2(0)p3(Ar) (1 + X^f^{e)) 

CRB (r) 

' 7r2rasNR'^'cos4(0)p3(iV) 

CRBS,(r) 



7r2ri)sNR'^'cos4(0)p3(iV) (1 + A2/2(e)) 
where fi{e) = ^^Jl^lf^^^^^^^^^ and f2{0) = ^2j2]^^lle)(N^-i) - P™"^ expressions, one can see that the two 
CRBs are quite similar for sources located in the central direction (i.e., for small 9 values) while at lateral directions 
(i.e., \9\ close to ^) the variable gain CRB is much lower than the equal gain one. This translates the fact that 
when the source location information contained in the time delay profile becomes weak, it is somehow partially 
compensated by the location information contained in the received power profile especially for the estimation of 
the range value. This observation is illustrated by the simulation experiments (cf. Fig. [5j Fig. [6] and Fig. |7]l given 
in Section IVll 

V. Near field localization region (NFLR) 

The radiating near field or Fresnel region is the region between the near and far fields lITTl corresponding to the 
space region defined by equation dSj. The latter depends on the source-antenna range, the signal wavelength, and 
the antenna aperture. 

Note that the lower bound of Fresnel region is related to the fact that in the immediate vicinity of the antenna, 
the fields are predominately reactive fields meaning that the E and H fields are orthogonal ifTTl . Therefore the 
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— ^— — ^— 1 should be kept out of the locaUzation region. However, the upper 
bound of the Fresnel region as given in equation ^ does not take into consideration the localization performance 
limit. For this reason, we suggest to define the near field localization region (NFLR) based on a target estimation 
performance relative to the application at hand. More precisely, if Stdmax > is the maximum standard deviation 
of the localization error that is tolerated by the considered application, i.e.. 



IE 



d-d 



< Stdn 



(34) 



where d = {x, y)'^ (resp. d) refers to the location vector (resp. its estimate), then the NFL region can be defined as 
the one for which the minimum standard deviation (given by the square root of the CRB) satisfies condition ([34]). 



Since E 
expressed as 



d-d 



E((x — x)^) + E((^ — y)^), the previous condition on the minimum MSE can be 



a/CRB(x) +CRB(y) < Stdr, 



(35) 



Now, the source coordinates can be rewritten according to 



X = r sm{9) = gx{r,6) 



y = rcos{6) = gy{r,0) 



and hence, by using the delta method in Ell , one can express 



CRB(x) + CRB(y) = Vglir,d)CVg.,{r,e) + Vg^{r,e)CVgyir,d) 



(36) 



where 



Vgx(r,0) 



dgx dgx 

dr do 

dPy dgy 

dr do 



[cos{0) rsvcL{0)Y 



[s\.u{9) rcos{0)y 



CRB (6') CRB(r,6 
CRB(6',r) CRB(r) 



A straightforward derivation of (1361 ) leads to 



CRB(x) + CRB(y) = r'^CRR{0) + CRB(r) 



(37) 
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and therefore the NFL region is defined as the one corresponding to 



Vr2CRB(6') +CRB(r) < Stdr, 



(38) 



An alternative approach would be to use a maximum tolerance value on the relative location error, i.e., a threshold 
value e such that 



\ 



E 



d-d 



< e 



(39) 



which corresponds to 



CRB{e) + 



CRB(r) 



< e 



(40) 



For example, in the conditional case, equation ( |40l ) becomes 

1 



2TD 



SNR 



GN{r,d)<e' 



(41) 



where 



GN{r,e) 



E,g{e)+E,g{r)/r^ 



From a practical point of view, equation (1411 ) can be used to tune the system parameters in order to achieve a 
desired localization performance. Different scenarii can be considered, according to the parameter, we can (or wish 
to) tune. 



Scenario 1: One can define the minimum observation time to achieve a desired localization performance at a given 
location and a given SNR value as 



GAr(r,( 



(42) 



SNR 



Similarly, one can also define the minimum SNR value for a target localization quality as 

Gjv(r,0) 



^SNR_„(^'' 



2e2r 



(43) 



In Section |Vll we provide simulation examples to illustrate the variation of these two parameters w.r.t. the source 
location. 
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Scenario 2: The previous parameters can be also be defined for a desired localization region (instead of a single 
location point (r, 6)) as 

Tmin{Rd) = max Tmin{r,9) (44) 
^SNR^J^d) = max Z)sNR_(r,0) (45) 
For example, if we are interested into the surveillance of a space sector limited by rmin < r < Vmax and 

— Gmax < 9 < Oraax, then 

" max Gat (r, 6*) 



I I VLI li 1-1 L 

262DSNR 

~ ^ o GN{r max, Gmax) (46) 

2e'^SNR 

^SNR„„„ = 

^ o orn^ ^iXmax 1 G-ijiax) (47) 

where the second equality holds from the observation that, away from the origin. Gat is a decreasing function w.r.t. 
the angular and range parameter^- 



Scenario 3: One can also wish to optimize the number of sensors with respect to a desired localization region Rd 
and for a target localization quality e. In that case, the minimum number of sensors needed to achieve the target 
quality can be calculated as 

Nmin = arg inin 

VI. Simulation 

In this Section, three experimental sets are considered. The first one is to compare the provided exact CRB 

expressions of lemma [T] with the CRB expressions given in |[T5l . Also, in that first experiment, we used the 

maximum likelihood approach (with the exact model in (dJ and ([2])) to validate our exact CRB derivation and to 

illustrate the gain in location estimation accuracy when using Q instead of (0]) in the data model. In the second 
''This is not an exact and proven statement but just an approximation tliat expresses the fact that the localization accuracy decreases when 
the source moves away from the antenna or towards its lateral directions. 
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experiment, we investigate the effect of considering the variable gain model instead of the constant gain model for 
near field source localization. Finally, the third experiment is to illustrate the usefulness of the near field localization 
region as compared to the standard Fresnel region. 

In all our simulations, we consider a uniform linear antenna with = 15 sensors and iter-element spacing 
d = ^ {\ = 0.5m) receiving signals from one near field source located at (r, 9). The sample size is T = 90 (unless 
stated otherwise) and the observed signal is corrupted by a white Gaussian circular noise of variance a^. In the 
conditional case, the source signal is of unit amplitude (i.e., a{t) = 1 Vt). 

The dotted vertical plots in all figures represent the upper and lower range limits of the Fresnel region given by 
©. 

A. Experiment 1: Comparison with existing work 

In Fig. |2j we compare the three CRB expressions for the source location parameter estimates versus the range 
values in the interval [0 50m] and versus ang le valueCl in the interval [-90° 90°]. The noise level is set to 
cr^ = 0.001 which corresponds to SNR = 30 dB. A similar comparison leading to similar results is given in Fig. 
[3] for a noise level set equal to cr^ = 0.5 which corresponds to SNR 3 dB. The source angle is ^ = 45° in the 
comparison versus range values and r = 20A in the comparison versus angle values. 

From these figures, one can observe a non negligeable difference between the exact CRB and the proposed one 
in |[T5l especially at low range values: i.e., the given CRB in |[T5l can be up to 30 times larger than the exact CRB. 
Also, from Fig. [21(c) and Fig. [3](c), contrary to the given CRB in |[T5l . the exact one varies with the range value 
with a relative difference varying from approximately 60% for small ranges to when r goes to infinity. 

From Fig. |2](b) - Fig. [3](b), one can observe that the lowest CRB is obtained in the central directions, this 

observation can be seen from the TE given in lemma |2] where the factor — \a\ is minimum for this directions and 

■— ' cos(y) 

goes to infinity when |^| |. 

Note that the provided Taylor expansion of the exact CRB is more accurate than the one obtained by expanding 

the time delay expression before CRB derivation i.e., the one in |[T5l . 

^Note that the curves w.r.t. the angle parameter (i.e., Fig. [2] Fig. [3] Fig.|7] Fig.[TO]and Fig. ll lb are not symmetrical around ^ = because 
we have chosen the first sensor for the time reference as shown in Fig. [T] 
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(a): e = 45° = 0.001 N = 15 T = 90 (b): r = 20X = 0.001 N = 15 T = 90 
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Range r Angle B 



(c): e = 45° 0^ = 0.001 N = 15 T = 90 (d): r = 201 = 0.001 N = 15 T = 90 



* Exact CRB in (9) 
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Fig. 2. CRB comparison: Exact conditional CRB versus approximate CRB in (15J in high SNR case 



(a): e = 45» = 0.5 N = 15 T = 90 (b): r = 201 <J^ = 0.5 N = 15 T = 90 
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(c): e = 45" 0^ = 0.5 N = 15 T = 90 ^ (d): r = 201 = 0.5 N = 15 T = 90 
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Fig. 3. CRB comparison: Exact conditional CRB versus approximate CRB in (15| in low SNR case 

To validate our exact CRB expression (see Fig. we have considered the Maximum likelihood (ML) algorithm. 
The ML function is computed for the equal gain conditional model, and its maximization is ensured by using 
Newton-Raphson method (6000 Monte Carlo runs are considered for the ML estimator MSB plot in Fig. IH. One 
can see that these CRB expressions can be reached and that we can gain up to 60% of MSB reduction by avoiding 
the standard model approximation in (01). 
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(a): = 45° 0^ = 0.001 N = 15 T = 90 
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Fig. 4. Exact conditional CRB validation 

B. Experiment 2: EG versus VG cases 

To this end, we have to ensure first that the received power in the two cases (i.e., constant and variable gain 
cases) is the same for the reference sensor (i.e., dividing the power of constant gain case per the square of the 
range as explained in Section IIV-AI ). To better compare CRB expressions, we consider two contexts where ^ = 0° 
for the first context (central direction) and = 85° for the second context (lateral direction). One can observe from 
Fig. |5]and Fig. [6] that, for small \0\ values, the constant gain CRB is quite similar to the variable gain one while 
at lateral direction (i.e., 16*1 close to ^) the variable gain CRB is much lower than the equal gain one, due to the 
extra information brought by the considered gain profile. 

This can be seen again from Fig. [T] where we can observe the large CRB difference for high |^| values. 

C. Experiment 3: Near Field Localization Region 

The plots in Fig. [8] represent the upper limit of the NFL region for different tolerance values. From this figure, 
one can observe that the Fresnel region is not appropriate to characterize the localization performance. Indeed, 
depending on the target quality, one can have space locations (i.e., sub-regions) in the Fresnel region that are out 
of the NFLR. Inversely, we have space locations not part of the Fresnel region that are attainable, i.e., they belong 
to the NFLR. 
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(a): 6 = 0° a'' = 0.001 N = 15 T = 90 



(b): 9 = 85- 0^ = 0.001 N = 15 T = 90 




* CRB with equal Gain 
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Fig. 5. CRB comparison: Equal Gain versus Variable Gain cases for range estimation 



(a): e = 0» 0'' = 0.001 N = 15 T = 90 



(b): 8 = 85° 0^ = 0.001 N = 15 T = 90 
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Fig. 6. CRB comparison: Equal Gain versus Variable Gain cases for angle estimation 



Fig. |9] compares the NFL region in the variable gain and equal gain cases with SNR = 30 dB. One can observe 

that in the lateral directions the NFLR associated to the variable gain model is much larger than its counterpart 

associated to the standard equal gain model. Also, in the short observation time context (i.e., Fig. |9](a)) the NFLR 

is included in the Fresnel region while for large observation time (i.e., Fig. |9](b)) the NFLR region is much more 

expanded and contains most of the Fresnel region^. 

^Except for the extreme lateral directions where the target quality can never be met since the CRB goes to infinity for \6\ — >■ ^. 
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(a): x = 2QX a" = 0.001 N = 15 T = 90 



(b): r = 20X = 0.001 N = 1 5 T = 90 
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Fig. 7. CRB comparison of the Equal Gain and Variable Gain cases versus angle value 
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Fig. 8. Near field localization regions for different values of the target quality: Stdmax 



In Fig. [To]- Fig. [TT] we illustrate the variation of the two parameters Tmin and ^g^R w.r.t. the source location 
parameters and for a relative tolerance error equal to e = 10%. From these figures, one can observe that -DsNR 
and Tjnin increase significantly for sources that are located far from the antenna or in the lateral directions. 



VII. Conclusion 

In this paper, three important results are proposed, discussed, and assessed through simulation experiments: (i) 
Exact conditional and unconditional CRB derivation for near field source locaUzation and its development in non 
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Fig. 9. Comparison of NFL regions in equal and variable gain cases 




<1 range 


= 10 


m 


* range 


= 25 


m 


range 


= 45 


m 




-100 -50 50 

Angle 8 



Fig. 10. Variation of the minimum observation time versus the source location parameters 



matrix form. The latter reveals interesting features and interpretations not shown by the CRB given in the literature 
based on an approximate model (i.e., approximate time delay), (ii) This CRB derivation is generalized to the variable 
gain case where the exact expression of the received power profile is taken into account. This generalization allows 
us to investigate the importance of the power profile information in 'adverse' localization contexts (i.e., for lateral 
lookup directions). 

{in) Based on the previous CRB derivation, a new concept of 'localization region' is introduced to better define the 
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E = 10% N = 15, T = 3000 



6 = 10% N = 15, T = 3000 
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Fig. 11. Variation of the minimum deterministic SNR versus the source location parameters 

space region where the locaUzation quaUty can meet a target value or otherwise to better tune the system parameters 
to achieve the target locaUzation quality for a given location region. 

Appendix A 
Proof of Lemma [5] 



A direct calculation of matrix Q in (|2TI ) using equation ^ leads to 

Q 



( \ 
Qi Q2 



where 



Qi 



Q2 



Q3 



( 

fee fer 

^ fr-e frr 

( , \ 



J 



J 



V2 V2 

v-sdiag{a a) 
Ot 



Or 

v-sIt 



the entries of Qi are given by (l22l ). (l23l) and (l24l ). (g) represents the Kronecker product, vi = ^(70 7)^7-51, 



V2 = ^{iQlfTr, 



^2 = ^7^7r.> 
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and V3 = II7II . 

Because the CRB of the range and the angle parameters is equal to the 2x2 top left sub-matrix of the inverse 



/ 



matrix Q , Schur lemma 11191 can be used and the results will be as Q 



\ 



where Qc = Qi 



X X 



After a straightforward computation, one obtain 

/ 

Qc = 



fe 



for 



fn 



'-iviV2 + v[v'^) frr - ^-^{V^ + V'i) J 



Now, by comparing this expression of Qc to the expressions in (l28l) . (|291 ) and (I30l ). one can rewrite 



leading finally to 



/ 



Qc - 2TD^^Yi 



\ 



CRBS,(r) 



E^gjr) 
det{Q^) 
1 



Evg{0) 



CRBU0) 



CRBl (r,9) 



2™SNr7 E,g{9)E,g{r) - E,g{r,ef 

det(Qc) 

1 \ E„g{r) 



2TD^^^J E,g{e)E,g{r) - E,g{r,e)^ 

Eygir,e) 

det(Qc) 

1 \ E,g{r,9) 



2™SNR/ E,g{6)E,g{r) - E,g{r,ef 
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Appendix B 



Proof of Lemma [6] 



For the unconditional case, the considered unknown parameter vector is ^ = (0, r, (T^ , u^)^ which leads to the 

/ \ 
Fi F2 



following 4x4 Fisher Information matrix FIM 



where the 2x2 matrices Fj are given by 



/ \ 

fee fer 



fr-e frr 



I \ 

feai fea^ 



J 



( \ 



By using Schur's lemma for matrix inversion |fT9l . one can obtain FIM ^ 



H 



where L = Fi 



FsFg^Fi^ 



u X 



. F3 and L are 2 X 2 matrices and their inverse can be computed easily as 



det 



u —X 



-X V 



( 



\ 



CRB(r) CRB(6',r 
\^ CRB(r,e) CRB(e) j 



where 



fr 

fe 

fri 



1 



deti 
1 

det I 
1 



(/6»(72Ci + fea^C2) 



det\ 



deti 

fa^cr^faie " fcr^a^fa^e 



(48) 
(49) 
(50) 
(51) 
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Cl = fa^a^-fa^-e-fa^^a^fa'^e (52) 
C2 = fa^^crlfa-^e-fcT-^a^Juie (53) 
C3 = fala^Ja^^r — fala^fa^r (54) 
C4 = /o-2o.2/o-2r — /o-2 0-2/0-2 r (55) 

det = uv — (56) 

Now, it remains only to compute the entries of the FIM by using equation (fTSl ) and taking into account that 
the matrix XI = agb{d,r)b{9,r)^ + cj^Iat and its inverse is given as S^^ = ^(Ia^ — ^b{9,r)b{6,r)^) where 

^= SNR + 11^11' ^"'^ = [7o,7ie^'"S--- ,77V-ie-''^"-r. 

A straightforward (but cumbersome) computation leads to 

2T 

foe = ^(l-SNR||7f)(7^7,)2 

-(l + SNRl|7||')((7 7)^-re)2 

+C SNR hf (||7,f + (70 7)^(f, re)) 

= _(i_SNR|l7f)(7^7.)' 

-(l + SNR||7||2)((7 7)^-r,)2 

+C SNR hf (||7,!|2 + (70 7)^(t, Tr)) 



Ire 



^(l-SNR||7f)(7^7.)(7^7.) 



-(1 + SNR |l7ir)((7 7)' re)((7 7)' t,) 

i>R 1171 

Tl 



-C SNR (7^ 7. + (7 if (re r,)) 

||4 



Cr4(C SNR)2 



W = ^(A^C2-||7ir(2C7-||7f)) 



/o-2(t2 



Cr4(C SNR)2 



= a^C^SNR^^ 
2T 

^"^^ ^ a^C^SNR^^'^^^) 

/ 2 = ^^(7-7 ) 
a2C2SNR^^ ^"-^ 

f = 2^ I T- ^ 



28 



By replacing these entries in equations (|48])-(|56]). we obtain 



u = 



V = 



X = 



2rSNR2 1 




|2 


(1 + SNR 


IItI 




2TSNR2 1 


|7| 


|2 


(1 + SNR 


IItI 
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2TSNR2 1 




|2 



(1 + SNR1I7I 



leading finally to Lemma [6] result 
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